The formation, ripening and stability of epitaxially strained island arrays 
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We study the formation and evolution of coherent islands on lattice mismatched epitaxially strained 
films. Faceted islands form in films with aniostropic surface tension. Under annealing, these islands 
ripen until a stable array is formed, with an island density which increases with film thickness. Un- 
der deposition, an island shape transition occurs, which leads to a bimodal island size distribution. 
In films with isotropic surface tension we observe continual ripening of islands above a certain film 
thickness. A stable wavy morphology is found in thinner films. 
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Coherent (dislocation-free) islands form to relieve the 
strain associated with lattice mismatched heteroepitaxial 
thin films. Their subsequent self-assembly into periodic 
arrays is of great interest as the arrays can be used to cre- 
ate quantum dot structures of importance in semiconduc- 
tor and optoelectronic devices. Such island arrays must 
have a narrow size distribution in order to be of use in 
applications. Of particular interest is whether the island 
arrays that form are energetically stable or metastable 
configurations that will ripen. Here we show in annealing 
simulations, that anisotropy in surface tension is neces- 
sary for the formation of a stable (roughly periodic) array 
with a narrow size distribution. Moreover, we show that 
the presence of a cusp in the surface energy is essential for 
reproducing the experimentally observed increase in is- 
land density with increasing film thickness. We also show 
that a single cusp in the surface energy (along with elas- 
tic relaxation) is sufficient in order to explain the island 
shape transition [1], which occurs in growth experiments, 
and the associated bimodal island size distribution. 

We study the evolution of an elastically isotropic sys- 
tem using continuum theory. The surface of the solid is 
at y = h(x, i) and the film is in the y > region with the 
film-substrate interface at y = 0. The system is modeled 
to be invariant in the z-direction, and all quantities are 
calculated for a section of unit width in that direction. 
This is consistent with plane strain where the solid ex- 
tends infinitely in the z-direction and hence all strains in 
this direction vanish. 

We assume that surface diffusion is the dominant mass 
transport mechanism, leading to the following evolution 
equation [2]: 
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where D s is the surface diffusion coefficient, r\ is the num- 
ber of atoms per unit area on the solid surface, f2 is the 
atomic volume, T is the temperature, kg is the Boltz- 
mann constant, s is the arc length and /x is the chemical 
potential at the surface. 

In our previous work [3,4] we showed that \x can be 
expressed as 
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where n is surface curvature, is the angle between the 
normal to the surface and the y-direction and j(6) — 
j(6) + d 2 j/80 2 is the surface stiffness (with ^(0) being 
the surface tension). Sijki are the compliance coefficients 
of the material, <7y is the total stress in the material, <7y 
is the mismatch stress in the zero strain reference state 
and fj,i\h) is the reference state free energy per unit 
length in the x-direction. The reference state is defined 
as a fiat film of thickness h confined to have the lateral 
lattice constants of the substrate. 

Linear stability analysis predicts that a flat film thin- 
ner than the linear wetting layer thickness, h c , is stable 
at all perturbation wavelengths and is marginally stable 
to perturbations of wavelength A c for thickness h c . The 
expressions for h c and A c are given in [3,4]. Above h c 
the flat film is unstable to a larger and larger range of 
wavelengths A_ < A < A + until for infinitely thick films 
the film is unstable to all perturbations of wavelengths 
larger than A = A c /2. 

We simulated the surface evolution given by Eqs. (1) 
and (2) using the numerical scheme described in our ear- 
lier work [4] . We used the cusped form of surface tension 
given by Bonzel and Preuss [5], which shows faceting in 
a free crystal at 0°, ±45° and ±90°. df^\h)/dh was ob- 
tained from ab-initio quantum mechanical calculations of 
Sii_ x Ge x grown on Si(00I) (for details see [6]). All our 
simulations start from a randomly perturbed flat film 
with an initial thickness denoted by C . 

When perturbations larger than a critical amplitude 
[3,4,6] are applied to a flat film, faceted islands develop in 
the film during both annealing and growth, as illustrated 
in Fig. 1. The film first becomes unstable at wavelength 
A ~ 507(0)/Me 2 , where M is the plain strain modu- 
lus and e is the lattice mismatch. The islands which 
form from this perturbation typically have a width of 
about 10% of the unstable wavelength. Both the critical 
wavelength and the faceted island widths scale as e -2 , 
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as observed in experiments [7-10] in which islands de- 
velop from long ripple like structures (corresponding to 
our model of plane strain). 

All results discussed henceforth refer to Ge/Si(001) 
though the same trends were seen in Geo.5Sio.5/Si(001). 
Islands form in a 'chain-reaction ripple' effect (i.e., is- 
lands tend to develop near other islands) as is illustrated 
in Fig. 1 . This mode of growth has also been observed in 
experiment [11,12]. The ripple effect occurs because the 
growth of the island destabilizes the flat film at its bound- 
aries. After initial island formation we observe island 
ripening occurring over much longer time scales (about 
50 times longer). 
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FIG. 1. Evolution of a random perturbation on a 25 mono- 
layer thick Ge film on a Si(001) substrate. In the first graph 
the dashed line is the initial perturbation, the thin solid line is 
the surface at t=0.005s and the thick solid line is the surface 
at t=0.032s. The second to sixth graphs show the surface at 
times t=0.044s, 0.068s, 0.094s, 0.123s and 2.181s. The final 
graph is the stable steady state island array. Note the ripple 
effect in island formation and the later island ripening leading 
to a stable island array. 



their sides are steeper than 45°. This shape transition is 
observed experimentally [1,14,17,18], and is sometimes 
referred to as the "pyramid-to-dome transition". The 
driving force behind it is the increased elastic relaxation 
experienced by tall-narrow islands. Theoretical equilib- 
rium calculations with isotropic surface tension [19] show 
a continuous increase in island aspect ratio with increas- 
ing island volume, as elastic effects dominate surface ten- 
sion effects. The sharp rather than smooth transition in 
growth mode we observe is due to the anisotropic nature 
of the surface tension and in particular the presence of a 
facet at 45°. Note that contrary to existing explanations 
of the island shape transition (see e.g. [20]), the transi- 
tion occurs without an additional facet orientation at a 
larger angle. 
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During annealing the islands are fully faceted. Their 
tops are faceted at 0° and their sides at 45°. This shape 
is preserved as the islands grow, i.e., the islands main- 
tain a fixed diameter-height ratio (as seen in experiment 
[7,13-15] and theory [16]). During deposition, on the 
other hand, an interesting transition is observed in the 
island shape. Initially, the islands are fully faceted as 
during annealing. However, when the islands reach a 
certain diameter, they stop growing laterally and only 
vertical growth occurs. This critical diameter is about 
40nm for Ge islands grown on a Si(001) substrate (during 
annealing we never observed islands which exceeded this 
diameter). Thus the islands become tall and narrow, and 



FIG. 2. Distribution of island cross-sectional area, A (re- 
call that islands are infinitely long in the z-direction), during 
directed deposition of Ge on a Si(001) substrate. The rate of 
deposition is 5.2nm/s, and the initial film height is 10 mono- 
layers. The dashed vertical line shows the separation between 
the early growth mode in which the island height-width ratio 
is preserved and the later vertical growth mode. 

The transition in island shape and growth mode is 
clearly reflected in the size distribution shown in Fig. 
2. Narrow island size and spacing distributions are seen 
during early deposition (see Fig. 2, 20 equivalent mono- 
layers). These narrow distributions are observed in many 
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experiments [1,7,14,15,18,21-23]. During later deposi- 
tion (30 equivalent monolayers) a bimodal distribution 
forms as some of the islands pass from the fully faceted 
to the tall-narrow shape. At later times (e.g. 50 equiv- 
alent monolayers) nearly all islands have the tall-narrow 
shape. At this stage the distribution becomes quite sym- 
metric and evolves at a fixed distribution width (increas- 
ing its mean). Similar results were observed in experi- 
ment [1,14,17,18]. 
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FIG. 3. Island density in a stable array of Ge islands on 
a Si(001) substrate after ripening has ended. C is the initial 
flat film thickness, and h m i is the thickness of one monolayer. 
The error bars refer to the span of island densities observed 
with different initial surface morphologies. 

One of our central observations is that annealing of a 
perturbed flat film with anisotropic surface tension leads 
to the formation of a stable array of islands. This result is 
consistent with several experimental systems [21,24-26], 
and is in contrast with films of isotropic surface tension 
where the islands ripen indefinitely. Theoretical studies 
also predict stable island arrays [27-30]. The crucial term 
in determining the stability of an island array apart from 
anisotropic surface tension and a film-substrate interac- 
tion is the presence of an elastic contribution due to is- 
land edges. This contribution is automatically present in 
our calculations and does not need to be introduced sep- 
arately. Theoretical works that ignore this term [22,31] 
predict continuous ripening. 

Our simulations show that the density of islands in the 
stable array increases with increasing film thickness (see 
Fig. 3). An increase in island density with film thickness 
has also been seen in many experiments [15,23,26,32-35]. 
Indeed Miller et al. [26] and Kamins et al. [35] performed 
annealing experiments and Leonard et al. [32] performed 
experiments with very small deposition rates. These 
three experiments clearly show the increase in island den- 
sity as film thickness increases. This result was predicted 
by Daruka and Barabasi [28] in minimal energy equilib- 



rium calculations. Here we show for the first time that 
the increase in island density also results from evolution 
simulations. This observation is particularly important, 
since other evolution studies [29] predicted a decrease in 
island density with increasing film thickness. We believe 
this is due to the smooth form of surface tension used 
in [29]. Indeed when we carried out simulations with a 
smooth form of surface tension similar to that used in 
[29], we also observed a decrease in island density. This 
clearly demonstrates the importance of using a cusped 
form of surface tension to accurately model evolution of 
faceting surfaces. 
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FIG. 4. Average island cross-sectional area, A (recall that 
islands are infinitely long in the z-direction) , in a stable array 
of Ge islands on a Si(001) substrate after ripening has ended. 
C is the initial flat film thickness, and h m i is the thickness of 
one monolayer. The error bars refer to the standard deviation 
in island sizes observed throughout all samples of the same 
film thickness. 

As can be seen in Fig. 4, the island size also shows 
a slight increase with increasing film thickness, with is- 
lands increasing in width from 25nm to 40nm, and cross- 
sectional area from lOOnm 2 to 500nm 2 (recall that is- 
lands are infinitely long in the ^-direction) . The island 
size at 7 monolayers is larger than expected due to fi- 
nite size effects. Note that even the smallest islands 
have a finite non-zero size. Experiments indeed see is- 
lands forming only above a certain size which increases 
with increasing film thickness [15,21-23,35]. However, as 
the experiments which were performed for annealing and 
showed stable island arrays tended not to vary the film 
thickness, it is difficult to compare our results with ex- 
perimental observations. Our result is in accordance with 
that predicted in equilibrium calculations by Daruka and 
Barabasi [28]. 

When surface tension is isotropic, corresponding to 
films above the roughening transition temperature, flat 
film evolution during annealing is very different from that 
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described above. Perturbations in films thinner than h c 
decay, and flat films with thickness h c < h < h c + A, 
where Awl monolayer, develop a stable smooth wavy 
morphology at A c . That is, perturbations of other wave- 
lengths decay and perturbations of wavelength A c grow to 
a finite amplitude. This is a mode of growth neither seen 
nor predicted before. Stable, non flat morphology, has 
previously only been predicted for faceting films [27-30]. 
In fact, other groups maintain that isotropic films should 
be unstable to ripening [29,30,36,37]. While films are lin- 
early unstable to perturbations of wavelengths A_ < A < 
A + , our simulations show that the nonlinearity stabilizes 
the growth of wavelengths close to A_ and A + . As a re- 
sult the growth of the perturbation saturates and stops 
at a finite amplitude, as seen by Spencer and Meiron [38] 
for infinitely thick films. When films are sufficiently close 
to the linear wetting layer thickness, the range of nonlin- 
ear saturation extends over the entire range of linearly 
unstable wavelengths and so a stable wavy morphology 
is observed. For films thicker than h c + A, initially a 
wavy structure at the most unstable wavelength forms. 
The hills of these waves then ripen on larger and larger 
length scales, until isolated islands are left that continue 
ripening. 
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